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We investigate asymptotic properties of long polymers grafted to convex cylindrical and spherical 
surfaces, and, in particular, distribution of chain free ends. The parabolic potential profile, predicted 
for fiat and concave brushes, fails in convex brushes, and chain free ends span only a finite fraction 
of the brush thickness. In this paper, we extend the self-consistent model developed by Ball, Marko, 
Milner and Witten to determine the size of the exclusion zone, i.e. size of the region of the brush 
free from chain ends. We show that in the limit of strong stretching, the brush can be described by 
an alternative system of integral equations. This system can be solved exactly in the limit of weakly 
curved brushes, and numerically for the intermediate to strong curvatures. We find that going from 
melt state to theta solvent and then to marginal solvent decreases relative size of the exclusion zone. 
These relative differences grow exponentially as the curvature decreases to zero. 



I. INTRODUCTION 

Long polymer chains grafted on a surface exhibit dif- 
ferent behavior from their random-walk state. Limited 
space available to the attached chains leads to their 
strong extension. The first model of grafted chains was 
pioneered by Alexander and de Gennes 0, 0, , who 
assumed that all chain ends are localized on the outside 
of the brush. Later Semenov p| and Milner, Witten, 
and Cates 0, 0, @ released the constraint on chain ends 
and developed a self-consistent field (SCF) theory of the 
brush. Similar a ppr oach was pursued by Zhulina and 

co-workers SEaia. The main conclusion of the SCF 

theory for flat brushes was that stretching of the chains 
is determined by a parabolic mean-field potential, and 
free chain ends span the whole area of the brush. 

Similar results hold for concave brushes. However, in 
convex brushes, direct minimization of brush free energy 
in the SCF approach leads to unphysical state with nega- 
tive density of chains. To explain this effect, Semenov Q 
has proposed existence of an "exclusion zone" , i.e. region 
of a brush free of chain ends. Later Ball, Marko, Milner 
and Witten have developed an integral equation ap- 
proach to a brush with exclusion zone that gave an exact 
solution for cylindrical brushes in a melt state. Subse- 
quently Li and Witten [13j have found exact solution for 
a cylindrical brush in marginal solvent in the limit of an 
infinite curvature. 

The size of the exclusion zone in spherical brushes, as 
well as in cylindrical brushes in solvents at intermediate 
curvatures is still c omp aratively ill-understood. Several 
simulations [ill liH ITa liH fl8| , semi-analytical calcula- 
tions [2, Il9j. and numerical approaches 0, |2(J have 
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been introduced. However, none of these fully exploits 
the limit of strong stretching to provide an analytical re- 
sult at intermediate curvatures, or gives a quantitative 
answer in the limit of very strong curvatures. 

In the present paper we investigate both cylindrical 
and spherical brushes in the melt, theta and marginal sol- 
vent states under the strong stretching assumption. We 
show that in this approximation brush properties may be 
described by a system of two integral equations. Com- 
pared with the results of Ball et all 0] , proposed equa- 
tions are linear and do not require self-consistent tuning 
of additional variables. In the limit of weak curvatures, 
these equation can be solved exactly for both melt and 
solvent states. We find that in this limit exclusion zone 
height exponentially decreases as radius of curvature in- 
creases. However, the rate of exponential decrease varies 
between solvent states, as well as between solvent and 
melt states. This suggests that the nature of the ex- 
clusion zone decrease is qualitatively different in various 
systems. 

In the region of intermediate to strong curvatures, the 
brush equations can be solved numerically. We find that 
in the limit of strong curvature exclusion zone occupies 
between a minimum of 42.4% of brush height for a cylin- 
drical brush in a marginal solvent to a maximum of 75.9% 
for a spherical brush in a melt. Finally, at all curvatures 
the relative size of the exclusion zone is larger in spherical 
than in cylindrical brushes. And in each of these systems 
exclusion zone decreases as one goes from melt to theta 
to marginal solvent. 



II. BRUSHES UNDER STRONG STRETCHING 

In the gaussian chain approximation, extension of a 
chain of N elements to length r costs ^(r 2 /Na 2 )(kT) in 
free energy. Here a is the Kuhn's length of the chain 
segment, and T is temperature. In the discussion that 
follows we will assume kT = 1. 

In grafted brushes, chains are additionally stretched 
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In a mechanically stable brush, osmotic pressure is 
pushing chain segments away from the grafting surface, 
so that n(x) is a monotonically decreasing function of x. 
Therefore, /x can be used in place of a: as a new coordinate 
(Fig. Pi. 

Following eq. (J2J , linear density of chain segments at 
point fi for the chains ending at point /iq is given by 



FIG. 1: Geometry of the brush, xq and Q are coordinate and 
potential at the exclusion zone boundary. S(x) is the area 
spanned by the chains at distance x away from the grafting 
surface. 
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by the repulsive forces from the adjacent chains, as 
monomers get expelled from the overfilled areas near the 
substrate. In the framework of a self-consistent field ap- 
proach, effect of adjacent chains can be described by an 
effective potential n(x). Physically this potential is re- 
lated to the osmotic pressure acting at a given point in 
space and is determined by an average concentration of 
chain segments at that point. A given /j(x) thus deter- 
mines the stretching profiles of the polymers and thence 
the monomer concentration profile. This profile must in 
turn be consistent with the original (i(x). Determining 
li(x) is thus a problem of self-consistency. 

The free energy of one chain in potential field /j,(x) is 
given by 
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The equilibrium properties of the brush are then deter- 
mined by the partition function Z — J^exp(— _F[a;(n)]), 
where summation is carried over all possible chain confor- 
mations. For brushes with high grafting density, chains 
are strongly stretched so that the partition function Z 
and other equilibrium properties are dominated by con- 
formations with lowest free energy Hence, in the 
strongly-stretched limit all chain are assumed to be at 
their lowest energy state. 

It was noticed by Semenov p| and Milner et al 0, Q, 
that the free energy has striking resemblance to the 
action of a classical particle moving in a potential field 
—fi(x). Hence the minimization problem is equivalent to 
finding a classical path of a particle moving in a given 
potential, with conservation of energy yielding 
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where constant of integration was chosen so that there 
is no pulling force on the open end of the chain: 



(dx/dn)\ 3 
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Now that we found chain extension in a given potential 
fj,(x), we can proceed to the self-consistent constraints im- 
posed on the potential. The first of these constraints is 
imposed by the geometrical restrictions on segment den- 
sity. In melt, for instance, density fluctuations are pro- 
hibitively expensive, and local monomer concentration 
is assumed constant throughout the brush. In solvent, 
on the other hand, the concentration may vary, but it 
is directly related to the chemical potential via the free 
energy of mixing |2l| . 

We start by introducing a function S(x) related to 
the geometry of the brush. It describes the area avail- 
able to chains at distance x from the grafting surface 
(Fig. [J), and, therefore, is presumed known. The phys- 
ically most interesting geometries can be summarized 
by S(x) = Cd+i{R + x) d , with d being "dimension" 
of the grafting surface (d = 0,1,2 for flat, cylindri- 
cal, and spherical surfaces respectively), and R being 
the radius of curvature of the surface. The constant 
Cd = d ■ 7r d / 2 /r(l + d/2) is the area of a unit sphere 
in d dimensions. Assuming that local volume fraction of 
monomers is 4>{x), the total area occupied by monomers 
at distance x is 4>(x)S(x). In the melt state, the volume 
fraction is constant <j){x) = 1, while in solvent state, </> 
varies with the distance. In the present paper, we will 
look at two types of solvents: a marginal solvent where 
(j, oc (j), and theta solvent where \i cx <p 2 21]. Consider- 
ation of other solvents, such as general case of good sol- 
vents, would require detailed analysis of local monomer 
interactions in and lies beyond the scope of this pa- 
per. Altogether <p(x) = r y s fi s (x), where s equals to 0, 1/2, 
and 1 for melt, theta and marginal solvents respectively, 
and 7 S is an unimportant proportionality coefficient. 

The total area <j>(x)S(x) spanned by the chains at dis- 
tance x from the grafting surface is comprised of the 
segments of the chains ending beyond x. Let a(x') be 
the total number of chains with free ends lying beyond 
x'. Then <t>{x)S{x) = a\ £ da{x r ) ^\ x ,^ x - Here m is a 
microscopic length close to segment packing parameter. 
Using fj, as a coordinate, the constraint on the segment 
density becomes 
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where S(fi) = S(x(fi)). The second line represents a well- 
known Abel equation and can be readily inverted 
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The second self-consistent constraint on the potential 
fi(x) is imposed by the length of the chains. The degree 
of polymerization for chains extended through point fi is 
given by 
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where P is potential at the grafting surface (Fig. QJ. 
Then condition of monodispersity would imply that 
N(fi) — No for all fi where chain ends are present. 

In a simpler case of concave and flat brushes, the two 
integral equations and © can be readily solved. In- 
deed, in those brushes chain ends span the whole vol- 
ume of the brush and N(fi) = N for all fi 6 [0, P). 
Then equation JfjJ immediately yields parabolic profile 
x oc (P — /i) 1 / 2 , so that chain distribution can be read- 
ily evaluated from (J5J. However, in the case of convex 
brushes, this straightforward approach results in physi- 
cally meaningless negative density of ends in the vicinity 
of the grafting surface. To avoid this unphysical behavior, 
one needs to introduce an "exclusion zone", i.e. a region 
of the brush near the grafting surface with no free ends 
inside. We denote the potential at the boundary of this 
exclusion zone as Q. Then the constraint on monodis- 
persity of chains reduces to N(fi) = N n for < fi < Q, 
while constraint on the absence of chain ends becomes 
<j{ijl) = (To for Q < fi < P. 

Introduction of these constraints significantly compli- 
cates the problem, as two equations © and J^Jl become 
interdependent. Still, in the melt state of cylindrical 
brushes, the two equation can be combined into single 
linear integral equation, that can later be solved exactly 
using electrostatic analogy 12]. In the presence of sol- 
vent or spherical brushes, on the other hand, the final 
equation becomes more complicated, and even non-linear 
in the case of a spherical brush, with no exact solution 
known to the authors. Additional complications arise 
from the singular behavior of the unknown dN/dfi and 
do/ dfi. 



To avoid this problem, we will rewrite equations (0 
and © in terms of one unknown function. The function 
S(x) describing available volume as a function of distance 
from the surface is presumed to be known. The function 
S(fi), on the other hand, is not. Hence, if we reformu- 
late and solve above equations in terms of this unknown 
function S(fi), then all the information about the system 
can be extracted from the comparison between S(x) and 
S(fi) = S(x(fi)). _ 

We start with introduction of the dimensionless poten- 
tials fi = fi/Q and P = P/Q, and proceed to rewriting 
equations an d © m terms of S(fl). Surface density 
of chains cr(/i) is constant throughout the exclusion zone 
(j2 > 1), so that equation Q can be rewritten as 
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where in the last step we substituted da /dfi' from © and 
integrated out one of the integrals. The new equation 
allows to find S(fi > 1) if S(fi < 1) is known. Immediate 
advantage of this new integral equation is its implicit 
incorporation of the exclusion zone and independence of 
Q or any other parameter of the system. 

Similar procedure can be applied to (jHJ. In this case 
N(fi) is constant outside the exclusion zone (fi < 1), and 

x(fx) = C d 1 ^ d S 1 / d (fi) — R. Once again we can combine 
© and its inverse so that 



S 1/d (fi < 1) = yjl-fi + a tan- 1 



S'/^fi') 




(8) 



[fi'-l] 1 ' 2 (fi'-fi) 



where a = [R 2 /2N 2 Qa 2 ] 1/2 = ^S^ d (P), and we have 
introduced rescaled S(fi) = (n 2 /8N§Qa 2 ) d ^ 2 S(fi). Obvi- 
ously, this equation allows to find S(fi < 1) if S(fi > 1) is 
known. Due to linearity, similar rescaling of S(fi) to S(fi) 
can be performed in equation Q , leading to a closed set 
of two integral equations Q and JHJ . Each one of these 
is linear, as opposed to the non-linear system of Refill. 
Additionally, the system implicitly incorporates ex- 
clusion zone so that the structure of the final solutions 
is not hidden by additional constraints (such as require- 
ment of positivity of end density [13]). 

The only parameter entering J7J) and JSJ is the scaled 
potential P which is related to the curvature of the graft- 
ing surface. In the melt state, this relation is immediately 
apparent from the conservation of brush volume. The 
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d+l 
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while the volume of a flat brush with the same number of 
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chains and the same grafting density is (d+ X)Cd+iR ho- 
Here ho is equilibrium thickness of the flat brush. Equat- 
ing the two volumes we get 
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In the presence of solvent, volume of the brush is no 
longer conserved and we need to use earlier definition 
of S(x) = C d+1 {ir 2 /8N 2 Qa 2 ) d / 2 {R + x) d . Since grafting 
surface (x = 0) corresponds to scaled potential jl = P, 
and open side of the brush (x = h) corresponds to jl = 0, 
we may write 
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refers to the flat brush, 
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1, as it follows from the equation JSJ in 



the limit of a flat brush (P — » +1). Hence R/h$ — 
(Q/Q//at) 1/2 5 1/d (P). The ratio Q/Q flat of the exclu- 
sion zone potentials can be found from the equality of 
grafting densities in curved and flat brushes. This den- 
sity is given by JSJ) 



S(P) 
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A. Weak Curvature Limit 

We now consider a limit of small curvatures. In this 
limit, exclusion zone height goes to zero, so that P — 
1 <C 1. Then function S{p) becomes nearly constant 
throughout the brush, and can be expanded about its 
value at the grafting surface Sp — S(P): S(jl > 1) = 
S P + AS(fi) and S(jl < 1) = S P + AS(/2). Finally, it is 
convenient to introduce new variable u = (1 — jl) 1 / 2 , for 
jl < 1, and u = [(fi - 1)/(P - 1)] 1/2 for jl > 1. 

To illustrate the approach, we start with the cylindrical 
melt brush (s = and d = 1). Then equations JJJ) and 
||SJ| transform into 
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Here e p — (P — l) 1 / 2 . The last two equations are linear, 
and can be combined. Keeping only leading terms in e p 
we get 
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The curvature of the brush is related to the quantity 
e p . Hence we are interested in the scaling behavior of Sp 
and A5' > (u) with e p . Applying standard separation of 
variables technique, we can regroup terms in (|15|l based 
on their u dependence and write down two separate equa- 
tions for Sp and AS>{u): 



In a flat brush P -> 1 and S(jl)/S(P) = (R + 
x(p.))/R — ► 1. Hence the radius of curvature and po- 
tential P are related via 



ln(e p ) 
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The remaining physical quantities can be immediately 
extracted from S(jl). The exclusion zone xq, i.e. the 
region of the brush near the grafting surface free of chain 
ends, can be found from (R + xq) 4 /(R + h) d — S(p, = 
T)/S(p, = 0). Hence 
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Here w is a constant which, in general, may depend on 
e p . The value of this constant is determined by the way 
original function S(u) is split between Sp and AS>(u). 
By construction, we want AS > (u) to be exactly zero at 
u = 1. Therefore substitution AS>(u) = — G(u) e p 
reduces Eq. (fT7|l to 



G(u) = -urn- u\n(u) + — G(v!) 



uln(^) du 1 
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with additional constraint that G(l) = 0. This equation 
has only one solution that satisfies given constraint. Us- 
ing results of Ref. 0, desired solution can be written in 
terms of the Meijer G-function pll2^|: 



G(u) = ~Gl 



3/2 3/2 
1 1 



(19) 



The constant ui introduced in l|16l) - (|17|) i s therefore also 
independent of e p and can be expressed in terms of the 
same Meijer G-function. Its numerical value is approxi- 
mately u> w 0.3864. Hence the final solution for a cylin- 
drical brush in a melt state is 



Sp = 



In. 



and 



A5> (u) = - G{u) e p . 

TT 
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As expected, when the brush becomes increasingly flat, 
Sp — > oo, and AS > — > 0. Hence, the radius of curvature, 
given by Eqs. © and (jllfl . and exclusion zone, Eq. (|12fl . 
become 
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(22) 



x Q /h » A5> (0) - A5>(1) = -ep. (23) 

These are the desired results for the cylindrical brush in 
a melt state. 

We now turn to the general case of a brush in a melt or 
solvent state. Once again, S(fi) can be expanded about 
its value near the grafting surface Sp, so that to the 
leading order in e p equations Q and JHJ reduce to 
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and 
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By analogy with the melt case discussed above, we sub- 
stitute (|2"5|) into l|2*4*jl and break up the resulting equation 
into two separate equations for Sp and AS>(m): 



AS>(u) 



2d ~ 



e p u [u> + In (it)] 



TT 2 , 



AS>(« 



; ,uln(5) 



(27) 



where iJ, is harmonic number 1251 . These solve to 
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Here G(u) and w are the same variables as those intro- 
duced in equations H18|) and (|19l) . Once again, as the 
brush becomes increasingly flat, Sp — > oo, and A5> — » 0. 
The radius of curvature (|llfl and exclusion zone l|12|) of 
a weakly curved brush are 
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Finally, it is convenient to write down a relation be- 
tween exclusion zone height and brush curvature. Us- 
ing earlier convention where d = 1 corresponds to the 
cylindrical brush, and d = 2 corresponds to the spherical 
brush, the exclusion zone height in the weak curvature 
limit is given by 



x Q /h ; 
xq/h : 
x Q /h : 



i e (fi>-i)-i*/fc» 



for melt 



for theta solvent 



for marginal solvent 



(32) 
(33) 
(34) 



where harmonic number H\i<x ~ 0.6137 and constant u> 
0.3864. 



B. Intermediate and Strong Curvatures 



1 gi/rf lXl + s) 1/2 1 
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For the general case of intermediate curvatures, we 
solve equations Q and © numerically. The unknown 
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cylinder, melt 
sphere, melt 
cylinder, e-solvent 
sphere, o-solvent 
cylinder, marginal solvent 
sphere, marginal solvent 




TABLE I: Relative size of the exclusion zone (xQ/h) at the 
limit of strong curvature for different brushes. 



- cylinder, melt 
sphere, melt 
cylinder, e-solvent 
sphere, B-solvent 
cylinder, marginal solvent 
sphere, marginal solvent 



R/h 



FIG. 2: Exclusion zone heights for cylindrical and spheri- 
cal brushes in melt and solvent states for a) strongly curved 
brushes, and b) weakly curved brushes. 



function S(p,) is proportional to the area spanned by- 
chains away from the grafting surface, and therefore 
is continuous and bounded everywhere on its domain. 
Hence we can approximate it by a piece-wise constant 
function and repeatedly iterate equations andJSJl un- 
til a stable solution is achieved. It is known |23 that 
iteration procedure for linear Fredholm equations with 
L,2 kernels converges when eigenvalue is small compared 
to the norm of kernel. Furthermore, Fredholm equations 
are also solvable when kernels have finite number of sim- 
ple poles. While these may not be used as a proof of 
convergences of the non-linear system (0 [SJ , the conver- 
gence may reasonably be expected. To further justify 
it, we have tried various initial states to verify that the 
system always iterates to the same final state. 

In the present work the function domain was split into 
4,000 segments, and iterations were repeated until ac- 
curacy of 10~ 9 was achieved. Use of piece-wise constant 
functions made it possible to analytically integrate out all 
the singularities, so that no further treatment of singular 
points was necessary. The results of numerical iterations 
are summarized in Figures |3 





Cylinder 


Sphere 


melt 


0.6340 


0.7594 


theta solvent 


0.4999 


0.6621 


marginal solvent 


0.4242 


0.5872 



TABLE II: Coefficients in the exponential dependence 
Ae~ aR/h ° of the exclusion zone height on radius of curvature. 
Predicted values are based on the weak curvature expansion, 
Eqs. 13211 - I34L and best-ht values are from numerical cal- 
culations, Fig. [5] 



Best-fit 
A a 



Predicted 
A a 



cylindrical brush, melt state 0.63 0.99 0.57 I 

cylindrical brush in theta solvent 0.52 1.56 0.49 tt/2 

cylindrical brush in marginal solvent 0.46 2.00 0.44 2 

spherical brush, melt state 0.71 0.52 0.57 1/2 

spherical brush in theta solvent 0.60 0.81 0.49 7r/4 

spherical brush in marginal solvent 0.52 1 .01 0.44 1 



III. DISCUSSION 

We have shown that exclusion zone is characteristic 
property of a brush in all considered solvent and melt 
states, for both cylindrical and spherical brushes. The 
size of an exclusion zone monotonically increases with 
curvature, until it reaches its finite maximum value at the 
extreme curvature limit. Our results are in perfect agree- 
ment with the known exact results for certain cylindrical 
brushes j 1 2L Il3j . Namely, we observe exact match with 
the exclusion zone height in the cylindrical melt calcu- 
lated for all curvatures by Ball et al Also, exclusion 
zone height matches with the exact result for a strongly 
curved cylindrical brush in marginal solvent calculated 
by Li and Witten However, it should come as no 

surprise, as all of these methods are based on similar SCF 
approach. 

Exact results for the weak curvature limit are summa- 
rized by equations - J51J. As one may expect, at 
weak curvatures dimension of the brush (i.e cylinder vs. 
sphere) is insignificant and exclusion zone height varies 
solely with the mean radius of curvature R/d. However, 
the strength of curvature dependence changes with sol- 
vent properties: decrease of the exclusion zone with ra- 
dius is fastest for marginal solvent, and is slowest for the 
melt. This might explain why it was not observed in ear- 
lier solvent simulations by Murat and Grest [17|, even 
though later simulations by Grest Q did indicate exis- 
tence of finite exclusion zone. Besides falling off faster, we 
also predict exclusion zone at extreme curvatures to be 
smallest for marginal solvents, and largest for the melt, 
with theta solvent in between (Table P). Comparison be- 
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FIG. 3: Volume fraction of chain segments in the case of 
extreme curvature (R = 0). All values are normalized to 
one at the boundary of the exclusion zone. In the order of 
increasing slope, the curves correspond to a) cylindrical brush 
in a theta solvent, b) cylindrical brush in a marginal solvent, 
c) spherical brush in a theta solvent, and d) spherical brush 
in a marginal solvent. 



tween analytical results in the weak curvature limit, Eqs. 
(|32|l - (|34|l and results of numerical calculations are sum- 
marized in Table [H] These results indicate very good 
agreement in the slope of the exponential decay. How- 
ever, weak curvature limit, i.e. limit of R — > oo, gives 
somewhat underestimated values of the intercept point 
in the depicted range of curvatures. 

In a spherical brush under marginal solvent conditions, 
the maximum size of the exclusion zone is xq ss 0.5872 h. 
This is significantly lower than 94% estimated by Li and 
Witten [lij using variational approach. However, such a 
difference is not surprising, as these authors have found 
that exact solution for the cylindrical brush is also sig- 
nificantly lower than the value predicted from variational 
approach. Additionally, lower size of the exclusion zone 
predicted in this work is in line with later simulations by 
Grest Il4| and numerical calculations by Wijmans and 
Zhulina |l9j. and Dan and Tirrell [20j . 

We find that the size of exclusion zone is always a 
monotonically increasing function of curvature , sim ilar to 
the cylindrical brushes under melt conditions {12] . How- 
ever, this contradicts with numerical calculations per- 
formed by Dan and Tirrell, who used diffusion-equation 
approach to the self-consistent problem |20(. These au- 
thors observed a sharp decrease in exclusion zone size at 
very strong curvatures, and it was found to be completely 
independent of any parameters of the system. However, 
these authors also observed unexpected drop in monomer 
density near the grafting surface. Hence one may spec- 
ulate that reported decrease in exclusion zone height is 
not a property of the diffusion-equation approach, but is 
rather a consequence of diverging numerical calculations 
near singularity. 

Still, possibility of a decrease in the exclusion zone 
height at large curvatures remains open. Stretching of 



TABLE III: Power law exponents for the dependence of vol- 
ume fraction on distance from grafting surface 4>(. x ) K x ° ■ 
Predicted values are based on equation 13511 . 





Best-fit 


Predicted 


cylindrical brush in theta solvent 


-0.4932 


- 1/2 


cylindrical brush in marginal solvent 


-0.6660 


- 2/3 


spherical brush in theta solvent 


-0.9871 


- 1 


spherical brush in marginal solvent 


-1.332 


- 4/3 



the chains continuously decreases away from the grafting 
surface, and fluctuations should play larger role in the 
outer region of the brush. As this effect becomes even 
more pronounced in strongly curved brushes, it may af- 
fect exclusion zone height and cause its decrease at large 
curvatures. We will leave this question for future re- 
search. 

The volume fraction 4>(x) of chain segments in differ- 
ent brushes in the limit of extreme curvature is shown in 
Figure El Deep inside the exclusion zone (x <C xq ) the 
volume fraction follows a power law, with corresponding 
best-fit exponents shown in Table IhTTl These fitting ex- 
ponents are in excellent agreement with predictions of 
a largely simplified model of the brush, which assumes 
that all chain ends lie on the outside of the brush. That 
model, first pioneered by Alexander and de Gennes 
El and later extended by Daoud and Cotton j2(| 
for spherical brushes, and Birshtein and Zhulina [27| for 
cylindrical brushes, predicts that the volume fraction of 
chain segments should scale with distance from the graft- 
ing surface as 

<t>{x) oc {R + x)- 1 -^- 1 ^ 2 ^ (35) 

where d is dimensionality of the brush (d = 1 for cylin- 
drical brush, and d = 2 for spherical brush), and the 
Flory exponent v depends on the solvent: v = 3/5 for 
good solvent, and v = 1/2 for theta solvent [l^. I n t ne 
vicinity of the grafting surface, i.e. far away from the ex- 
clusion zone boundary, results of our model match with 
Eq. (|35|l . as shown in Table ITTT1 However, as one ap- 
proaches the boundary of the exclusion zone, deviations 
from scaling models arise, and the power law exponent 
becomes smaller than that predicted by l|35|l . Outside 
exclusion zone, the density profile falls off much faster 
and no longer follows the power law. 

Similarly, in the brushes with non-zero radius of curva- 
ture R, simplified model of Daoud and Cotton becomes 
inapplicable and deviations from density profile (|35|l in- 
crease. Qualitatively, absolute value of the power law 
exponent decreases as radius increases. 

Finally, the density profiles for the chain free ends are 
shown on Figure The presence of the exclusion zone 
does restrict chain ends to a fraction of a brush, but it 
does not qualitatively affect their distribution: at the 
boundary of the exclusion zone the density of ends con- 
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FIG. 4: Density of free ends in a spherical brush a) at strong 
curvature limit, and b) in marginal solvent at different curva- 
tures. 



ever, away from the grafting surface, as one approaches 
boundary of the exclusion zone, the deviations arise. 
And, ultimately, at the opposite side of the brush, i.e. 
near the open surface, the brush properties approach 
those predicted by Semenov Q and Milner et al |y, Ll t| • 

Analysis given in appendix shows that approximation 
of strong stretching is suitable for most brushes in the 
limit of large degree of polymerization. Hence it should 
provide a good qualitative guide to the real systems, and 
in particular to the colloidal particles or comb polymers 
in a solvent environment. Detailed knowledge of chain 
conformations in these brushes may provide new mecha- 
nisms of controlling properties of these systems. 

Finally, the approach developed in this paper may 
prove useful in analyzing brushes with inhomogeneous 
structures, and, in particular, brushes with polydisperse 
chains or multi-component systems. Investigation of 
chain conformations in these systems, as well as anal- 
ysis of fluctuations due to non-asymptotic stretching, are 
some of the possible areas of future research. 
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tinuously approaches zero, while at the open surface free 
ends distribution maintains the same structure as found 
in flat and concave brushes. 



APPENDIX: APPLICABILITY OF THE STRONG 
STRETCHING APPROXIMATION. 



IV. CONCLUSIONS 

In conclusion, we have shown that exclusion zone occu- 
pies a finite fraction of the brush in all convex geometries 
considered. In both cylindrical and spherical brushes, in 
either solvent or melt state, the size of the exclusion zone 
is a monotonic function of curvature, continuously in- 
creasing as radius of curvature goes to zero. In the limit 
of strong curvature it reaches the maximum value, which 
is always a finite fraction of the brush thickness (42% to 
76% of the brush, Tabled). 

At weak curvatures exclusion zone height decreases ex- 
ponentially with radius. However, rate of this decrease 
depends on the solvent /melt state of the system, indicat- 
ing qualitatively different nature of the underlying mech- 
anism. 

In the vicinity of the grafting surface the properties of 
the brush are very well described by the scaling model of 
Daoud and Cotton |2(| . The match becomes nearly per- 
fect in the limit of strong curvatures (Table ITTT|) . How- 



The approximation used in this paper was based on two 
assumptions. First of all, grafted chains were assumed to 
be significantly stretched compared to their free radius of 
gyration. This allowed us to neglect fluctuations in chain 
trajectories. On the other hand, chain extension was as- 
sumed to be significantly smaller than the full extended 
length to ensure validity of the gaussian approximation 
|JT|| . Introduction of these constraints imposes some lim- 
its on the applicability of the theory, which are discussed 
below. 

Our analysis here is similar to that of Ref. IT2I so we 
will briefly summarize the main conclusions and extend 
them to the broader class of solvent states considered in 
this paper. 

The length of a chain in a dense brush is of the same or- 
der as the brush thickness h. Therefore, aforementioned 
constraints on chain length can be written as 

aN 1 / 2 < h < aN, (A.l) 
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TABLE IV: Scaling of brush height with grafting density and 
chain polymerization for strongly curved brushes. 





cylinder 


sphere 


melt 


(T l/2 jV l/2 


CT l/3jVl/3 


theta solvent 


^1/3^2/3 


a l/4 jV l/2 


marginal solvent 


(T l/4 jV 3/4 


a l/5 jV 3/5 



where N is chain polymerization, and a is characteristic 
size of one monomer. Here we are interested in asymp- 
totic behavior at large N. 

Let us first consider flat and weakly curved brushes. 
In a melt state flat brush has a thickness h oc Na 3 a/A, 
where a is the total number of grafted chains, and A is 
total area of the grafting surface. Therefore, constraints 
(lA.ljl are asymptotically satisfied by brushes with low 
grafting density. In the solvent state, brush height can 
be estimated from minimization of average free energy 
per chain F = h 2 /2Na 2 + w v (jfV/a, where v describes 
solvent state (y — 2 for marginal solvent, and v = 3 
for theta solvent), V — Ah is the total volume of the 
brush, and w v is unimportant proportionality coefficient 
related to the virial coefficients. The volume fraction of 
monomers in a flat brush is <f> = aa 3 N/V, and minimiza- 



tion of free energy gives h oc (cr/A)^ -1 " ( v+x > jV. Sim- 
ilarly to the melt case, nearly flat brushes in a solvent 
also asymptotically satisfy strong stretching approxima- 
tion at low grafting density. 

Similar considerations may be applied to the curved 
brushes. Considering the limit of strong curvatures, the 
scaling laws for the brush heights take the form shown in 
Table llVl Evidently, most brushes once again satisfy Eq. 
(IA.III : the cylindrical and spherical brushes in marginal 
solvent, as well as cylindrical brushes in theta solvent 
always obey (jA.ljl in the asymptotic limit of large N. The 
cylindrical brushes in a melt state, and spherical brushes 
in theta solvent asymptotically satisfy l|A.l|) . provided 
that their grafting densities are sufficiently large. Here 
large densities mean that the total number of chains a ^> 
Wg^a 11 / 8 for the sphere, and linear density of chains a 3> 
a -1 for the cylinder. 

The only system that violates (|A.I|I in the asymp- 
totic limit of large N is the spherical brush in a melt 
state. Specifically, this system fails the first of the two 
constraints, and therefore requires detailed accounting of 
different chain trajectories. Still, our theory should cor- 
rectly describe this brush in the limit of weak curvatures, 
and may give qualitative insights for the case of interme- 
diate to strong curvatures. 
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